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A  Look-Up  Perimeter  Intrusion  Sensor: 
Evaluation  of  a  Continuous  Array 


1.  INTRODUCTION 

This  report  describes  and  analyzes  a  look-up  sensor  intended  to  encircle 
the  perimeter  of  a  large  installation  and  detect  low  altitude  airborne  threats 
such  as  helicopters  or  hang  gliders. 

A  number  of  intrusion  sensors  have  been  developed  to  detect  intruders  on 

the  ground.  Their  operation  has  been  based  on  several  different  physical  phe- 

1-3 

nomena  including  acoustic,  seismic,  optical,  infrared,  and  radio  waves. 

Most  of  these  have  very  limited  detection  zones  or  require  the  intruder  to  make 

physical  contact  with  the  sensor.  In  contrast,  the  detection  fields  of  radar-like 

radio  frequency  systems  can  be  made  to  extend  significant  distances  from  the 

4 

sensor  and  can  thus  detect  an  intruder  where  conventional  sensors  can  not. 
(Received  for  publication  18  April  1984) 
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A  sketch  of  a  conceptual  look-up  sensor  is  depicted  in  Figure  1  where  the 
sensor  is  deployed  around  the  perimeter  of  an  installation.  The  detection  zone 
extends  upwards  from  the  ground  and  an  object  penetrating  this  zone  would  be 
detected.  However,  the  value  of  the  sensor  would  be  greatly  enhanced  if  it  could 
also  determine  the  position  of  the  intrusion  along  the  perimeter.  One  way  to 
accomplish  this  is  to  drive  the  sensor  with  a  pulse,  rather  than  a  continuous 
wave.  The  pulse  travels  along  the  sensor  and  excites  a  section  of  the  sensor 
proportional  to  the  duration  of  the  pulse,  thus  producing  an  upward  pointing 
beam.  The  beam  moves  along  the  sensor  and  intercepts  the  target.  Some  of 
the  energy  reflected  by  the  target  is  received  by  the  sensor  and  propagates  back 
to  its  input.  The  time  At  that  elapses  between  the  transmission  of  the  pulse  and 
the  return  of  its  reflection  from  the  target  is  a  measure  of  the  total  distance  to 
the  target.  The  elapsed  time  is  related  to  the  target's  position  by 

At  =  2{(X/v)  +  (Y/c)}  . 


where  X  is  the  distance  along  the  sensor  to  the  point  of  intrusion,  Y  is  the  alti¬ 
tude  of  the  target,  and  v  and  c  are  the  propagation  velocities  of  the  pulse  along 
the  sensor  and  in  free  space  respectively. 

The  generic  sensor  configuration  that  will  be  analyzed  consists  of  a  continu¬ 
ous  array  of  radiators,  each  coupled  to  a  suitable  transmission  line.  When  the 
pulse  propagates  along  the  transmission  line  the  group  of  elements  spanned  by 
the  pulse  forms  an  upward  pointing  beam  that  moves  along  the  sensor  at  its 
propagation  velocity. 

The  equations  describing  the  electromagnetic  properties  of  the  array  making 
up  the  sensor  structure  will  be  developed  by  applying  phased  array  theory.  How¬ 
ever,  since  the  targets  can  be  in  the  near-field  of  the  array  and  the  amplitude 
coefficients  of  successive  elements  can  be  different,  the  equations  that  result 
are  not  as  compact  as  the  far-field  approximations  often  associated  with  linear 
arrays.  The  derived  sensor  equations  will  then  be  used  as  the  basis  for  calcu¬ 
lating  the  radar  performance  of  the  system  from  which  some  preliminary  designs 
can  be  developed.  Finally,  the  effects  of  tapering  the  coupling  to  insure  a  con¬ 
stant  detection  sensitivity  and  the  temporal  shaping  of  the  pulse  to  control  the 
effective  sidelobes  of  the  beam  pattern  will  be  addressed. 
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2.  COMPUTATION  OF  SENSOR  FIELDS 
2.1  Dipole  Field* 

The  e  component  of  electric  field  intensity  due  to  a  thin  half -wavelength  long 
dipole  radiator  oriented  along  the  z  axis  (Figure  2)  can  be  written  as 

ikr 

g  *  V  — _  cos[(t  cos,e)/2]  x  (1) 

6  or  sin  6 

as  long  as  the  distance  r  from  the  center  of  the  dipole  to  the  field  point  exceeds 
a  few  wavelengths.  5  The  free  space  wavenumber  is  k  =  2*/ A,  Vq  is  the  complex 
amplitude,  9  is  the  polar  angle,  and  9  is  the  unit  vector. 


V 


Figure  2.  Dipole  Element  and  Its  Coordinate  System 


5.  Jasik,  H.  (Ed.)  (1961)  Antenna  Engineering  Handbook,  McGraw-Hill,  N.  Y. 
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The  computations  to  follow  may  be  further  simplified  if  9  is  set  equal  to  tt/2 
thus  restricting  attention  to  the  xy  plane.  In  this  case  Eq.  (1)  becomes 


(2) 


where  Q  is  seen  to  be  equal  to  -z. 

The  array  is  made  up  of  a  number  of  similar  dipoles  arranged  as  indicated 
in  Figure  3.  The  local  field  of  each  dipole  is  given  by  Eq.  (2)  in  terms  of  the 
coordinates  of  that  element.  In  contrast,  the  field  quantities  of  the  array  will 
usually  be  given  in  terms  of  the  coordinates  of  the  first  excited  element  or  the 
center  of  the  array.  All  dipoles  are  in  the  xz  plane,  lie  parallel  to  the  z  axis, 
and  are  separated  by  a  distance  d.  Each  element  is  considered  to  be  excited  by 
a  transmission  line  with  propagation  characteristics  described  by  an  attenuation 
factor  a  and  a  phase  constant  £i.  The  length  of  line  between  elements  is  d1,  that 
between  the  input  and  the  first  element  is  x^,  and  that  between  the  last  element 
and  the  matched  load  terminating  the  line  d^.  Primed  quantities  refer  to  the 
transmission  line  lengths  and  unprimed  quantities  to  the  spatial  locations  of  the 


Figure  3.  Sketch  of  Dipole  Array 


elements.  The  total  length  of  the  sensor  is  L  and  includes  the  length  x^.  The 
degree  of  coupling  (power)  between  the  transmission  line  and  the  nth  element  is 


The  physical  properties  of  the  array  need  not  be  specified  further  at  this 
time,  since  the  sensor  could  be  implemented  in  a  number  of  ways.  For  example, 
the  dipoles  could  be  proximity-coupled  to  a  long  transmission  line.  Instead,  the 
performance  of  the  generic  array  shown  in  Figure  3  will  be  investigated. 


2.2  Dipole  Array 

The  total  field  at  a  point  P(X,  Y),  due  to  the  linear  array  of  dipoles  depicted 
in  Figure  3,  is  the  sum  of  the  contributions  from  each  of  the  elements.  For  N 
sources,  this  sum  is 


E0(X,Y) 


N 

-E 

n=  1 


e[-j  (<ot-krn)]  „ 


(3) 


where  is  the  complex  amplitude  of  the  nth  element. 

The  power  density  is  found  from  the  average  Poynting  vector 
P  =  1/2  Re(E  X  H*)  to  be 


P  =  Re 


N  N 

EL 

n=l  m  =  l 


*  Jjk(rn-rm)]  .  . 

VnVm  2».r„r—  cx*n. 
o  n  m 


(4) 


where  zq  is  the  intrinsic  impedance  of  free  space.  The  normalized  field  intensity 

EoN  and  the  normalized  power  density  PN  are  obtained  by  dividing  Eg  and  P  by 

the  magnitude  of  their  respective  maximum  values. 

If  the  observation  point  is  sufficiently  far  from  the  array,  then  Eqs.  (3)  and 

(4)  can  be  specialized  and  some  additional  quantities  defined.  From  Figure  3  it 

is  seen  that  rn  in  the  numerator  of  Eq.  (1)  satisfies  the  relation  rR+^  =  rn  " 

(x  , ,  -  x  )  cos  <t>,  that  r,  =  R,  and  that  r  in  the  denominator  can  be  replaced 
n  *  l  n  l  *1  o 

by  R.  These  changes  are  valid  for  R  S  2[(N  -  l)d]  /X  where  (N  -  l)d  is  the 
length  of  the  array.  Dropping  the  vector  notation  and  setting  (*n+1  "  xn)  =  d- 
let  us  write  Eq.  (1)  as 


v  ■  'H,"*-kE’1  E 

n=l 


=l-j(n-l)kd  cos 


(5) 
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and  the  normalized  far-field  intensity  as 


E  =  J-j(ut-kR)]  y  e[-j(n-l)kd  cos  $]  /  y 
eFN  /  -»  n  /  /  i  n 


The  power  density  may  be  written  as 

PF  =  E0F  E0F^2zo  (7) 

and  the  normalized  power  density  as 

PFN  =  E0FN  E0FN  -  (8) 

If  V  =  V,  then  Eq.  (8)  reduces  to  the  familiar  expression  for  the  far-field 
power  pattern  of  a  linear  array  of  N  isotropic  elements  (or  dipole  elements,  as 

g 

in  this  case,  when  the  field  is  measured  only  in  the  xy  plane).  That  is. 


(9) 


The  normalized  power  patterns  PpN  for  several  arrays  with  element  spacing 
d  =  x/2  and  =  Vq  are  shown  in  Figure  4.  The  maximum  gain  occurs  at  <$>  =  90° 
(broadside)  and,  generally,  the  longer  the  array,  the  narrower  the  main  beam 
and  the  lower  the  average  sidelobe  level. 

The  ultimate  purpose  of  this  analysis  is  to  determine  the  radar  detection 
performance  of  the  sensor  as  a  function  of  the  target  coordinates  X  and  Y. 
Therefore,  it  is  more  appropriate  to  present  the  power  pattern  of  the  array  a.' 
a  function  of  the  X  coordinate  of  the  target  as  it  moves  along  the  length  of  the 
sensor  in  the  xy  plane  with  the  altitude  Y  taken  as  some  arbitrary  but  fixed 
value.  Placing  the  center  of  the  array  at  xq  =  L/2  yields 

2  2  1/2 

r  =  Hx  -  x  r  +  Y  ]  (10) 

O 


6.  Krauss,  J.D.  (1950)  Antennas,  McGraw-Hill,  N.  Y. 
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Figure  4.  Normalized  Far-Field  Power  Patterns  for  Several  Arrays, 
d  =  x/2.  (a)  N  =  16,  (b)  N  =  32,  and  (c)  N  =  64  (Contd) 


and 


$  =  tan_1[Y/(X  -  xQ)] 


(ID 


The  expression  for  the  normalized  far-field  intensity  obtained  by  inserting  Eqs. 
(10)  and  (11)  into  Eq.  (5)  and  then  normalizing  to  the  magnitude  of  the  maximum 
field  intensity  at  X  =  xQ  is 


E0FN 


2  2  \ 

ut-ki<x-x0r  +  y^]  I 


v  f.  (-j(n-l)kd  cos{tan_1[Y/(X-x0)]}  ) 
n  ' 

2  o  1^2 

Kx-x0)2  +  y2]  £  vn 

n=l 


(12) 
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(b) 


DISTANCE  ALONG  SENSOR  (m) 


Figure  5.  Normalized  Far -Field  Power  Density  at  a  Constant  Height 
Above  Sensor.  L  =  1  km,  x_  =  L/2,  d  =  X/2,  X  =  1  m,  V  =  V,  and 
Y  =  2(N  -  I)2d2/X.  (a)  N  =  32  and  (b)  N  =  64  (Contd) 


Re  (Efl  XH.) 

5  -  _ 0  <P 

*  '  2 


(15) 


When  Y  is  sufficiently  large,  Eq.  (15)  reduces  to  the  expression  given  by  Eq.  (7) 

since  then,  <t>  -  <t>  for  all  m.  .  /o 

m  2  2 

For  a  constant  height  target,  r  =  [(X  -  x  )  +  Y  )  , 

,  m  m 

<t>m  -  tan  1[Y/(X  -  x^)],  xm  =  x^  +  (m  -  l)d,  and  Xj  =  {L  -  (N  -  l)d}/2,  where 
L  is  the  total  length  of  the  sensor.  Similar  expressions  apply  for  the  index  n. 
With  these  substitutions,  the  power  density  expressed  by  Eq.  (15),  normalized 


to  its  value  broadside  to  the  array,  is  plotted  in  Figure  6  for  three  values  of  Y 
to  show  the  effect  of  the  target's  close  proximity  to  the  array  elements.  Com¬ 
parison  of  the  patterns  shows  a  filling  of  the  near-in  nulls  when  the  target  is  not 
in  the  far -field  of  the  array.  This  is  as  expected  because  of  the  quadratic  phase 


error  and  the  increasing  amplitude  error.  The  apparent  narrowing  of  the  main 
beam  results  from  the  fact  that  for  a  target  near  the  ground,  the  angle  0  changes 
very  rapidly  with  X  compared  to  a  higher  target.  Thus  the  effective  linear 
resolution  of  such  a  system  will  be  a  function  of  the  height  of  the  target  as  well 


(c)  DISTANCE  ALONG  SENSOR  (m) 


Figure  6.  Normalized  Near-Field  Power  Density  at  a  Constant  Height 
Above  Sensor.  L  =  1  km,  xQ  =  L/2,  N  =  32,  d  «  >/2,  X  =  1  m,  and 
Vn  =  V.  (a)  Y  =  2  (N  -  l)2d2/X,  (b)  Y  =  (N  -  l)2d2/X.  and 
(c)  Y  =  0.  2(N  -  l)2d2/X  (Contd) 

as  the  length  of  the  array.  The  pattern  in  Figure  6(a)  should  be  compared  with 
that  in  Figure  5(a)  to  show  that  Eqs.  (15)  and  (7)  produce  progressively  closer 
results  as  R,  the  distance  from  the  center  of  the  array  to  P,  increases. 

2.3  Element  Excitation 

The  field  distribution  expressed  by  Eq.  (3)  should  be  modified  to  account 
for  the  phase  and  amplitude  of  the  signal  exciting  successive  elements  of  the 
array.  This  can  be  done  by  letting  take  the  form 


V  =y  gf-o'+jPHxi  +  (n-l)d']  K(1_|K|2)[(n-l)/2] 
n  o 


(16) 


where  Vq  is  the  transmission  line  input  voltage,  a  and  (3  its  propagation  con¬ 
stants,  and  K  the  element  coupling  coefficient.  With  these  changes,  Eq.  (3) 
can  be  written  as 


[(^*+j0)xi  -  jut] 


*VVoe 


N 

Y  el(-a+j0)(n-l)d'  +  jkrn]  K( x_ | K 1 2)(n-l)/2 

£l  ^  (17) 


and  the  power  density  found  with  Eq.  (15).  The  expression  for  the  far-field 
intensity  corresponding  to  Eq.  (17)  is 

r  ((^^xl  -i(<ot-kR)1  V  e{l('°f+jl3)d,'jkdCOS^(n‘1)} 

^flF  ■  vo  e  Z-f  R 

n=l 


K(l-|K|2)(n'1)/2  .  (18) 

From  the  preceding  discussion  it  can  be  seen  that  the  power  radiated  from 
successive  elements  decreases  steadily  because  of  the  attenuation  of  the  trans¬ 
mission  line  and  the  energy  coupled  to  the  previous  elements.  To  provide  equal 
power  radiated  from  each  element,  the  coefficients  must  be  varied  along  the 
array.  The  new  coefficients  Kn  that  will  ensure  that  equal  power  is  coupled  to 
each  element  are  computed  as  follows. 

The  power  reaching  the  position  of  the  nth  element  along  the  transmission 
line  can  be  expressed  as 


for  n  >  1,  where  Zq  is  the  input  impedance  of  the  transmission  line.  The  product 

PnKnK*  thus  is  equal  to  the  power  radiated  from  the  nth  element  and  Pn(l  "KnK*) 

is  the  power  available  to  continue  down  the  line  to  excite  the  remainder  of  the 

elements.  The  requirement  for  an  equal  amplitude  signal  to  be  radiated  from 

each  element  then  gives  the  recursion  formula  P  K  K*  =  P  ,  .K  ,  ,K* , ,  from 
°  n  n  n  n+1  n+1  n+1 

which  an  expression  can  be  constructed  that  gives  the  value  of  the  coupling  coef¬ 
ficient  for  any  element  in  terms  of  the  first.  The  result  is 
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(20) 


for  n  >  1.  The  power  remaining  immediately  after  the  last  element  (n  *  N)  is 
=  Pjj(1  "  anc*  *s  absorbed  by  the  load  terminating  the  feeding  cable. 

Successive  application  of  Eq.  (19)  under  the  restriction  of  Eq.  (20)  yields 
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where  PT  =  VqVo/2Zo<  the  input  power.  Equation  (21)  is  solved  for  KjKj  and 
substituted  into  Eq.  (20)  to  obtain  the  remaining  Kn>  The  results  are  plotted  in 
Figure  7  for  several  values  of  the  parameter  A  *  pl/pt' 
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Figure  7. 
tation.  d' 


Element  Coupling  Coefficients  for  Uniform  Element  Exci- 
=  0,73  m,  A  =  1  m,  N  =  2000,  and  a  =  0.  0025  Neper/m 


19 


2.4  Pulse  Illumination 


When  a  pulse  propagates  down  the  transmission  line,  a  certain  time  must 
elapse  before  the  energy  reaches  points  farther  down  the  line.  Thus,  each 
element  of  the  array  is  energized  at  a  slightly  later  time  than  the  previous 
element  and  during  this  transient  interval  the  beam  forming  properties  of  the 
array  are  incomplete. 

To  analyze  this  behavior  it  is  assumed  that  the  transmission  line  and  the 
radiating  elements  are  not  dispersive  so  that  each  frequency  component  of  the 
pulse  is  similarly  affected.  To  account  for  the  propagation  time  along  the  trans¬ 
mission  line,  Eq.  (16)  is  written  as 


Vn'V>sf  ['-T-I  VoKne 


((-a+j|3Kxi  +  (n-l)d1)] 


n  <hkpi2> 


where 


f(t  -  x^/v)  =  1  for  0  <  t  -  xn/v  <  r 


=  0  otherwise. 


Thus  a  pulse  does  not  arrive  at  the  point  x'  on  the  transmission  line  for  x'/v  sec 
after  it  is  applied  to  the  input  and  persists  for  only  t  sec,  the  pulse  width,  after 
which  time  the  voltage  remains  zero  until  the  arrival  of  another  pulse.  The 
excitation  of  each  element  of  the  array  occurs  d'/v  sec  after  the  previous  element 


and  the  energy  radiated  by  each  element  reaches  the  field  point  P(X,  Y)  rn/c  sec 


after  that.  To  account  for  this,  Eq.  (17)  is  written  as 
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where 


f[t  -  (x^/v)  -  (rn/c)]  *  1  for  0  <  t  -  (x^/v)  '  (rn/c>  -  T 


0  otherwise. 


In  the  limit  as  t  -  f[t  -  (xn/n)  -  (rn/c))  —  1  and  Eq.  (23)  reduces  to  the  ex¬ 
pression  for  Eg  given  by  Eq.  (17).  In  a  similar  manner,  the  far-field  intensity  is 
represented  as 


E  -  V  et<-c*+j/3)xi  ‘  j(wt-kR>l 
0F  oe 
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17  - 6 -  •  (24) 

p=i 

The  last  two  equations  can  be  written  in  much  simpler  form  if  the  constraint 
that  equal  power  be  radiated  from  each  element  is  enforced.  Then  the  coupling 
coefficients  of  each  element  increase  to  compensate  exactly  for  radiation  and 
line  loss  and  Eqs.  (23)  and  (24)  can  be  written  as 
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Inspection  of  Eq.  (23)  shows  that  not  all  N  elements  of  the  array  may  be 
excited  simultaneously.  For  example,  when  vr  <  (N-l)d',  the  span  of  the  pulse 
is  less  than  the  length  of  the  array  and  the  coefficients  f[t  -  (x^/v)  -  (rn/c)]  will 
be  identically  zero  for  some  of  the  elements  and  the  maximum  gain  of  the  array 
will  not  be  realized.  Only  when  vr  ^  (N-l)d'  will  all  elements  be  illuminated  at 
the  same  time  for  at  least  part  of  the  pulse  duration  period.  When  the  array  is 
incompletely  illuminated,  the  maximum  number  of  excited  elements  is 

N'  =  VT/'d1  +  1. 

The  peak  value  of  Pp  =  (Re  Egp  X  H^p)/2  normalized  to  the  full  array  value 
is  plotted  in  Figure  8  to  show  the  effect  of  the  transient  phenomenon  on  the  power 
density.  Inspection  of  the  curves  shows  that  until  the  pulse  spans  the  entire  array 
the  beam  is  incompletely  formed  and  the  power  density  at  P(X,  Y)  is  below  its 
maximum  possible  value.  With  the  fill  time  F  *  (N-l)d'/v  and  r  =  0.  5F,  only 
half  the  array  elements  are  excited  and  the  greatest  gain  reached  is  6  dB  below 
the  full  array  gain.  For  r  =  F,  the  full  gain  is  obtained  only  for  the  instant  that 
the  pulse  remains  centered  on  the  array  and  every  element  is  excited.  For  r  >  F, 
the  gain  remains  at  its  maximum  value  for  proportionally  longer  intervals. 


Figure  8.  Power  Density  for  an  Incompletely  Excited  Array.  N  =  32, 
L  =  1000  m,  X  =  1  m,  v  =  0.73c,  d'  =  0.73  m,  d  =  X/2,  and 
F  =  (N  -  l)d'/v 
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2.5  Continuous  Array 

When  the  number  of  elements  on  the  transmission  line  is  large  enough  so 
that  vt  <  (N-l)d’,  the  array  is  continuous  and  the  width  of  the  pulse  on  the  line 
determines  the  length  of  the  excited  segment.  As  the  pulse  travels  down  the 
transmission  line,  the  position  of  this  subarray  also  moves  along  with  the  pulse. 
Thus  a  beam  is  formed  that  also  moves  down  the  line  at  a  velocity  equal  to  the 
propagation  velocity  of  the  line.  The  power  density  at  a  point  above  the  sensor 
is  a  time -dependent  function  that  grows  to  a  maximum  value  as  the  traveling 
beam  intercepts  it  and  then  decays  to  zero  gradually  after  it  has  passed. 

The  power  density  produced  by  a  pulse  traveling  down  the  continuous  array 
can  be  obtained  from  Eqs.  (25)  and  (26)  by  summing  over  every  element  of  the 
sensor  instead  of  the  few  that  made  up  the  finite  array.  To  demonstrate  the 
effect  of  this,  the  power  density  at  a  point  halfway  down  a  1-km  long  sensor 
populated  with  dipoles  along  its  entire  length  is  plotted  in  Figure  9  as  a  function 
of  time  for  several  field  point  heights.  The  abscissa  in  Figure  9  begins  at: 
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Figure  9.  Power  Density  for  Continuous  Array,  r  =  10.333  X  10*®  sec, 
L  =  1000  m,  A  =  1  m,  v  =  0.73c,  d'  =  0.  73  m,  d  =  \/2,  xi  =  0,  N'*32, 

X  *  L/2,  and  T  =  [xi  +  (N1  -  l)d'/2]/v  +  [Y2  +  (X  -iN'  -  l)d/2)2  ]1/2/c. 
(a)  Y  =  2(vT)2/ A,  (b)  Y  =  (v//A,  and  (c)  Y  =  20(vT)Z/A 
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the  time  required  for  the  center  of  the  excited  subarray  of  N'  elements  to  reach 
the  first  element  and  the  signal  transmitted  from  the  leading  element  to  have 
reached  PCX,  Y).  The  buildup  and  decay  of  the  density  traces  out  the  subarray's 
radiation  pattern.  However,  it  is  apparent  that  although  the  curve  has  many  of 
the  properties  of  the  patterns  shown  in  Figure  6,  its  sidelobe  structure  is 
markedly  different.  The  sidelobe  spacing  during  the  buildup  phase  is  smaller 
than  that  during  the  decay  phase  of  the  pattern.  This  is  easily  explained  and  is 
a  result  of  the  high  propagation  velocity  of  the  subarray  along  the  sensor.  It 
might  be  expected  that  the  contributions  at  P(X,  Y)  from  successive  elements 
would  arrive  at  time  intervals  of  d'/v.  However,  the  spatial  path  lengths  from 
successive  elements  become  progressively  shorter  up  to  the  element  immediately 
under  the  field  point,  after  which  they  become  progressively  longer.  Thus  during 
buildup,  contributions  from  successive  elements  reach  the  field  point  P(X,  Y)  at 
intervals  shorter  than  d'/v.  As  a  result,  the  elements  appear  to  be  more  closely 
spaced,  the  sidelobe  spacing  decreases,  and  the  sidelobe  decay  rate  increases. 
When  the  subarray  is  in  the  neighborhood  of  the  target  location  X,  this  effect  has 
disappeared  because  the  spatial  pathlengths  from  successive  elements  to  the 
target  are  nearly  equal.  In  this  region  the  pattern  structure  and  subarray  gain 
correspond  to  an  array  of  N'  =  vr/d'  +  1  elements.  Beyond  this  region,  the 
arrival  times  of  successive  contributions  are  delayed  because  the  spatial  path- 
lengths  are  progressively  longer.  The  result  is  that  the  subarray  appears  to 
have  fewer,  more  widely  spaced  elements.  Its  pattern  then  shows  greater  spac¬ 
ing  between  sidelobes  and  a  slower  decay  rate  than  during  the  buildup  phase. 

The  number  of  elements  excited  at  any  one  time  is  always  the  same  (once 

the  entire  pulse  is  on  the  sensor)  and  equal  to  N'  =  vr/d'  +  1.  However,  as 

discussed  above,  the  apparent  number  of  elements  observed  at  P(X,  Y)  varies  as 

the  subarray  propagates  across  the  sensor.  The  apparent  number  of  subarray 

elements  M  is  equal  to  the  number  of  element  locations  x  for  which 

n 

0  ^  t  -  (x^/v)  -  (rn/c)  £  t  .  (27) 

The  value  of  M  so  determined  is  plotted  on  the  right  in  Figure  9  for  each  of  the 
curves.  The  staircase  character  of  the  curve  is  a  result  of  the  integer  counting 
that  results  from  applying  Eq.  (27).  The  count  can  vary  by  one  element  for  the 
same  interval  depending  on  the  relative  placement  of  the  count  window.  It  is  this 


.v v'.v.vvyv .v.’-v 


'vyi'.v.v''. 


TI 


r-y,  vvr/\«  cj-  r Ji" JC"^>  P>*> 

ft 


variation  in  M  that  produces  the  noisy  structure  in  the  sidelobe  region  of  the 
pattern.  Notice  that  at  the  peak  of  the  beam  M  =  N'  =  (vr/d')  +  1  =  32  as  expected. 

The  effect  of  height  Y  above  the  sensor  on  the  power  density  is  apparent 
from  the  comparison  of  Figures  9(a),  9(b),  and  9(c),  which  should  also  be  com¬ 
pared  with  the  corresponding  curves  of  Figure  6.  The  structure  of  these  curves 
would  have  been  observed  at  P(X,  Y)  if  that  array  had  been  moved  slowly  along 
the  ground  over  the  distance  occupied  by  the  sensor. 

An  estimate  of  the  apparent  number  M  of  radiators  excited,  as  a  function  of 
time,  can  be  obtained  from  the  relation 
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M  =  {t/  (dt/ drn) }  +  1 


(28) 


whenever  N  »  N'.  The  expression  for  dt/dr  is 


dt  .  _d| _ (d/c)[X  -  (n-l)d] 

dr„ 


v  {Y2  +  [X-(n-l)d]2}1/>2 


(29) 


For  X  =  L/2  and  N  odd,  dt/drR  at  n  =  1,  n  =  (N  +  l)/2,  and  n  =  N  is  respectively 


dt  _  d^_  (d/c)(L/2) 
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d^  v 


(Yz  +  (L^/4)  I 


dt/dr(N+l)/2  =  d’/v  ’ 


(30) 


and 


dt  d1  ,  (d/c)(L/2) 

dr„,  v  Z  , 1/2 
[ Y^  +  (L  / 4) ] 


N 


For  reasonable  values  of  L  and  Y,  the  effect  ot  v  «  c  is  to  make  the  second  term 


in  Eq.  (30)  negligible  so  that  dt/drn  =  d’/v  for  all  n.  This  corresponds  to  moving 


a  physical  subarray  of  N1  elements  slowly  along  the  sensor  while  recording  the 
power  density  at  P.  To  show  this,  Eq.  (28)  is  plotted  in  F.gnre  10  for  several 
values  of  v. 

In  radar  calculations,  it  is  the  maximum  value  of  power  density  that  deter¬ 
mines  detection  performance.  Thus  only  the  peak  value  of  the  main  beam  need 
be  computed  to  find  the  gain  of  the  subarray.  This  occurs  when  the  subarray  is 
centered  under  the  target,  M  =  N’. 
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Figure  10.  Apparent  Number  of  Elements  Excited.  L  =  1000  m, 
X  =  1  m.  d'  =  0.73  m.  d  =  X/2,  and  Y  =  2(N'  -  l)2d2/x 


An  expression  for  the  gain  of  the  array  sensor  in  the  plane  0  =  x/2  can  be 
obtained  from  Eqs.  (25)  and  (26)  by  setting  fit  -  (x^/v)  -  (rn/c)j  =  1.  By  defini 
tion,  the  gain  G(0,<t>)  =  4xR2  P(e,tf)/P0.  where  Pq  is  the  same  input  power  that 
established  the  power  density  P (0,<t>)  and  R  the  distance  from  the  center  of  the 
array  to  the  field  point  P(X,  Y).  Setting  PQ  =  VoV*/2Zq  and  <t>  =  x/2,  we  obtain 
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where  rp  =  [Y2  +  (X  -  d{l(N'+l)/2]  -  p})2]1/2  and 

<t>  -  tan-1  [Y/(X  -  d{[(N'  +  l)/2]  -  p})]  with  a  similar  expression  for  the  index  q. 

p  a  2  i 

The  maximum  value  of  the  far-field  gain  is  =  4ir  N'  exp[-2axj  ]  for 

0d'  =  2ir.  By  reciprocity,  the  sensor  gain  on  receive  is  the  same  as  on  transmit. 


3.  RADAR  PERFORMANCE  ANALYSIS 


3.1  Computation  of  Received  Radar  Signal 


In  this  section,  the  performance  of  an  array  of  dipoles  used  as  a  radar 
sensor  will  be  analyzed.  The  system  consists  of  a  continuous  array  of  dipoles 
along  a  perimeter  of  length  L. 

The  power  appearing  at  the  terminals  of  the  transmission  line  due  to  a 
reflection  from  a  target  at  P(X,  Y)  is  computed  from  the  radar  range  equation 

7 

in  the  form 

PREC  =  {p0GT/4,rR2}  W4^2}  U2gr/M  <  (32) 


where  P  is  the  input  power;  G^,  and  GR  the  transmit  and  receive  antenna  gains 
respectively;  X  the  wavelength;  and  a  the  radar  cross-section,  assumed  here  to 
be  isotropic.  The  factors  in  the  first  bracket  give  the  power  density  at  the  target 
produced  by  an  antenna  of  gain  G^  and  input  power  PQ.  The  second  bracket,  when 
multiplied  by  the  first  gives  the  power  density  at  the  receiving  antenna  produced 
by  a  reflection  from  a  target  of  radar  cross  section  a  square  meters.  The  factors 
in  the  last  bracket  represent  the  effective  area  A  of  the  array  that  is  related  to 
its  gain  by  the  expression  G^  =  4irA/\  . 

Thus,  Eq.  (32)  may  be  used  to  compute  the  power  received  at  the  input  of 
the  system  from  the  interaction  of  the  transmitted  signals  with  a  target  entering 
the  detection  zone  of  the  sensor.  The  performance  of  the  sensor  system  depends 
on  the  magnitude  of  PRE(~,  relative  to  the  system  noise  as  will  be  discussed  in 
the  next  section. 

3.2  Probability  of  Detection 

The  quantity  that  defines  the  performance  of  any  radar  is  the  signal-to-noise 
ratio  S/N  at  its  detector.  For  this  analysis  S  =  and  the  noise  is  N  =  kTgB, 

where  k  is  the  Boltzmann  constant,  Tg  the  system  noise  temperature  in  °K,  and 
B  the  noise  bandwidth.  The  receiver  noise  at  the  detector  input  is  Rayleigh 
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distributed  and  the  received  signal  statistics  depend  on  the  size  and  shape  of  the 
target,  its  velocity,  the  operating  frequency,  and  the  pulse  repetition  frequency. 
In  the  present  case,  the  received  radar  signal  will  be  assumed  to  be  nonfluctuat¬ 
ing,  a  reasonable  assumption  for  target  sizes  on  the  order  of  a  wavelength. 

In  addition  to  the  internally  generated  and  galactic  antenna  noise,  most  radar 
receivers  must  also  contend  with  clutter  signals  from  the  ground  or  other  sources 
within  the  same  range -azimuth  cell  as  the  target.  A  variety  of  techniques  have 
been  developed  to  separate  the  desired  Doppler  shifted  signal  that  has  been  re¬ 
flected  by  the  moving  target  from  the  unshifted  signal  reflected  from  stationary 
clutter  sources.  As  a  result,  the  clutter  power  is  reduced  to  the  system  noise 
level  and  often  need  not  be  considered. 

The  signals  that  finally  reach  the  detector  are  the  receiver  noise  and  the 
target  signals,  the  clutter  signal  having  been  eliminated  by  appropriate  filters. 

The  probability  of  detection  PQ  and  the  probability  of  false  alarm  F^  depend  on 
the  threshold  at  which  the  detector  is  set,  the  S/N  ratio,  and  the  statistics  of 
the  target  signal.  The  that  can  be  expected  for  an  assumed  F^  can  be  cal¬ 
culated  using  well-known  methods  or  can  be  obtained  from  curves  that  have  been 

g 

prepared  for  a  large  range  of  system  parameters  and  fluctuation  statistics. 

For  convenience,  a  few  curves  appropriate  for  the  detection  of  a  nonfluctuating 
target  against  Rayleigh  noise  are  shown  in  Figure  11. 

3.3  Performance  Calculations 

At  this  point  it  is  necessary  to  select  a  set  of  system  parameters  upon  which 

to  base  the  performance  calculations.  The  pulse  repetition  frequency  (prf)  was 

selected  to  achieve  an  unambiguous  detection  of  a  target  at  1-km  altitude  along 

a  1  km  long  sensor.  Thus  prf  <  l/[(2Ld'/dv)  +  (2Y/c)]  <  50  kHz.  The  pulse 

width  was  chosen  to  provide  an  altitude  resolution  s  =  cr/2  of  50  m,  requiring  a 

pulse  width  r  55  0.  35  ps.  This  in  turn  gives  N'  =  vr/d'  +  1  »  77  for  the  number 

of  elements  simultaneously  excited.  To  improve  the  S/N  ratio,  100  pulses  are 

integrated.  Because  of  the  50-kHz  prf,  a  target  moving  at  the  moderately  high 

speed  (for  this  application)  of  20  m/sec  will  have  changed  position  by  0.  04  m 

over  the  100-pulse  presumming  interval.  Finally,  the  wavelength  is  chosen  to 

2 

be  \  =  1  m  and  the  radar  cross  section  to  be  0  dBsm  (1  m  ). 

With  the  noise  bandwidth  B  set  equal  to  1/r,  the  S/N  ratio  was  evaluated  as 
a  function  of  Y  for  several  values  of  P  .  The  PD  was  obtained  from  Figure  11 
and  the  results  plotted  in  Figure  12. 

Inspection  of  the  curves  show  that  a  peak  input  power  of  about  10  W  will 
provide  a  probability  of  detection  of  about  0.  9  along  the  entire  length  of  the 


8.  Skolnik,  M.I.  (Ed.)  (1970)  Radar  Handbook,  McGraw-Hill,  N.Y. 
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Figure  11.  Probability  of  Detection  as  a  Function  of  Signal-to-Noise  Ratio  for  a 
Nonfluctuating  Target.  Data  for  curves  obtained  from  Reference  8 


sensor  to  an  altitude  of  1  km.  The  short  interval  of  reduced  performance  at 
about  250  m  (for  PQ  =  1  W)  is  a  result  of  the  varying  near-field  gain  as  a  function 
of  altitude.  Over  this  narrow  interval,  the  decrease  in  S/N  produced  by  the  l/R4 
loss  is  greater  than  the  increase  in  the  square  of  the  near-field  gain.  The  rate 
of  increase  in  the  gain  decreases  with  altitude  until  it  reaches  the  far-field  limit 
of  Eq.  (31b). 

3.4  Measurement  of  Range  and  Altitude 

When  the  system  is  configured  in  the  monostatic  mode  as  shown  in  Fig¬ 
ure  13(a),  only  the  elapsed  time  between  the  transmission  of  a  pulse  and  the 
reception  of  its  echo  can  be  measured.  From  the  sketch,  it  is  clear  that  the 
elapsed  time  &t  is  related  to  the  target  location  by 
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TARGET  ALTITUDE  (m) 

Figure  12.  Probability  of  Detection  for  a  Target  at  Various  Altitudes. 
ar=  0.  0025  Neper/m,  Xj  =  0  m,  KjKj  =  -45  dB,  a  -  0  dBsm,  r  =  0.  35 
psec,  and  X  =  1  m 


At  =  2(Xd'/vd  +  Y/c)  (33) 

and  that  the  target  can  be  anywhere  along  the  straight  line  described  by 
Y  =  cXd'/vd  +  cAt/2.  Thus,  the  measurement  cannot  distinguish  between  range 
and  altitude  and  only  the  total  range -altitude  delay  can  be  measured. 

To  overcome  this  problem,  the  system  can  be  configured  into  the  bistatic 
mode  as  shown  in  Figure  12(b).  Here  another  receiver  Rx  2  has  been  added  to 
measure  the  difference  At'  between  the  transit  time  Ld'/vd  of  the  pulse  directly 
to  Rx  2  and  the  propagation  time  of  the  pulse  down  the  line  to  the  target,  up  to 
the  target,  and  back  down  to  the  line  to  Rx  2,  The  propagation  time  for  this  is 
Xd'/vd  +  2Y/c  +  (L  -  X)d'/vd.  Thus  At'  is 

At'  =  2Y/c  ,  (34) 

and  gives  for  the  target  altitude 


2  [( X  d 1  /  v  d )  +(Y/c)] 


target 


From  Eq.  (33)  the  range  is  found  to  be 

X  =  vd(M  -  At')/2d'  .  (36) 

The  addition  of  the  second  receiver  allows  the  unambiguous  measurement  of 
both  the  target's  range  and  its  altitude. 

It  is  interesting  to  use  Eq.  (33)  and  the  curve  in  Figure  9(a)  to  compute  the 
target  range  X.  From  Figure  9(a),  the  time  between  T  and  the  arrival  of  the 
pulse  at  P(X,  Y)  is  about  792.  5  d'/v.  To  this  must  be  added  that  portion  of  T 
that  accounts  for  the  propagation  time  of  the  contribution  from  r  16,  which  is 
tlg  =  {Y2  +  [X  -  (N'  -  l)d/2]2  }*/2/c.  Thus  ^t  is  twice  the  sum  of  792.  5  and  tjg. 
Taking  Y  =  2(N'  -  l)2d2/\  there  results  X  =  vd[(At/2)  -  (Y/c)]/d'  =  499.  9  *  500  m, 
the  known  range  of  the  target. 

3.5  Pulse  Shaping  for  Lower  Sidelobes 

All  patterns  shown  thus  far  were  calculated  assuming  that  the  excitation  for 
each  element  was  the  same.  Such  uniform  arrays  produce  the  maximum  gain 
obtainable  for  a  given  number  of  elements  and  exhibit  a  relatively  high  sidelobe 
level.  Figures  4  and  5  show  that  the  peak  of  the  first  sidelobe  is  about  13  dB 
below  the  peak  of  the  main  beam.  The  sidelobe  structure  of  an  array  is  greatly 
affected  by  the  relative  excitation  levels  of  its  elements.  The  average  sidelobe 
level  as  well  as  the  distribution  of  the  levels  of  the  sidelobes  can  be  controlled 
by  choosing  the  appropriate  excitation  distribution  (taper)  of  the  array.  Com¬ 
monly  used  distributions  are  the  Taylor  and  Chebyshev  tapers.  The  former 
produces  a  pattern  that  reduces  the  first  n  sidelobes  to  a  uniform  level  and  then 
allows  the  remainder  to  fall  off  gradually.  The  latter  produces  a  sidelobe  dis¬ 
tribution  with  all  sidelobes  at  the  same  reduced  level. 

In  conventional  arrays,  the  excitation  of  the  elements  is  controlled  by  the 
construction  of  the  beam  former,  which  is  designed  to  produce  the  desired 
taper.  In  the  present  case,  however,  the  excitation  levels  cannot  be  incorporated 
into  the  construction  of  the  array  because  the  subarray  is  formed  by  the  pulse 
as  it  moves  down  the  line.  Thus  it  is  necessary  to  control  the  shape  of  the  pulse 
by  suitably  modulating  its  amplitude.  To  demonstrate  the  principle,  a  simple 
linear  modulation  that  produces  a  triangular  pulse  on  a  pedestal  is  used  instead 
of  the  more  complicated  Taylor  or  Chebyshev  tapers.  The  time  modulation  will 
result  in  an  element  excitation  distribution  that  can  be  represented  by  the  expres¬ 
sion 


=  1  -  i #  m  =  lto 


M  (M  odd) 


vm  =  VM+l-n  =  1  '  (M~2MKl"*i~  *  ^  =  1  to  M/2  (M  even)  , 


where  u  is  the  fractional  height  of  the  pedestal. 

The  effect  of  modulating  the  excitation  pulse  can  be  computed  with  Eq.  (2  5) 

if  a  factor  v  is  inserted  to  control  the  amplitudes  of  the  element  contributions 
m 

for  which  f[t  -  (x^/v)  -  (rn/c)]  =  1.  With  this  change  Eq.  (25)  becomes 


E  =  V  K  eH<*+j0)xi  "  i*1*1 

e  o  l 


V  ft-' 

m  l 


Jj3d'(n-1)+ jkrn  /_0 


where  1  s  m  <  M  with  M  determined  from  Eq.  (27)  or  Eq.  (28).  m  =  1  is  made 
to  correspond  to  the  first  element  for  which  f[t  -  (x^/v)  -  (rn/c)]  2:  0.  The  pat¬ 
terns  that  result  when  v  from  Eqs.  (37)  or  (38)  is  substituted  into  Eq.  (39)  are 

m 

shown  in  Figure  14  for  three  values  of  pedestal  height.  It  is  apparent  that  the 
patterns  in  Figures  14(a)  and  14(b)  exhibit  a  lower  sidelobe  level  and  broader 
beamwidth  than  does  that  of  the  uniform  array  in  Figure  9(a). 

The  principal  benefit  of  lower  sidelobes  in  radar  is  to  reduce  the  possibility 
of  interfering  signals  entering  the  radar  receiver. 


4.  CONCLUDING  REMARKS 


A  continuous  array  of  radiators  appropriately  coupled  to  a  transmission  line 
is  proposed  for  a  look-up  sensor  to  detect  low -altitude  airborne  intruders.  The 
chief  advantage  of  such  a  system  is  that  it  produces  a  continuous  vertical  zone  of 
protection  that  extends  all  the  way  to  the  ground  with  no  dead  zones  or  shadows 
and  that  follows  the  perimeter  of  the  installation.  To  achieve  this  complete 
coverage  zone  however,  the  physical  sensor  must  also  encircle  the  entire  instal¬ 
lation. 

The  use  of  horizontal  dipoles  produces  maximum  gain  in  the  vertical  direc¬ 
tion.  The  detection  performance  decreases  gradually  [Eq.  (1)]  away  from  verti¬ 
cal  going  to  zero  toward  the  horizon(s)  away  from  the  longitudinal  vertical  plane. 
For  the  present  case,  the  detection  performance  is  symmetrical  about  the  longi¬ 
tudinal  vertical.  In  certain  applications,  it  may  be  desired  to  decrease  the 
detection  sensitivity  toward  the  inside  of  the  installation  and  to  increase  the 
sensitivity  toward  the  outside.  In  such  cases,  radiators  with  appropriately 
shaped  element  patterns  chosen  to  achieve  the  proper  coverage  could  be  used. 

In  fact,  the  element  patterns  of  the  radiators  could  be  varied  along  the  sensor 
to  shape  the  coverage  as  needed. 

The  maximum  length  of  the  sensor  is  primarily  limited  by  the  attenuation  of 
the  feeding  network.  To  extend  its  length,  line  amplifiers  could  be  used  to  raise 
the  signal  level  in  the  feed  transmission  line.  The  coupling  coefficients  of  the 
sensors  must  be  decreased  to  offset  the  gain  of  the  amplifier  in  order  to  maintain 
a  uniform  detection  zone.  Also,  the  allowable  pulse  repetition  frequency  will  be 
affected  by  the  length  of  the  sensor. 
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